home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
Loadstar 174
/
174.d81
/
star locator
(
.txt
)
< prev
next >
Wrap
Commodore BASIC
|
2022-08-26
|
4KB
|
149 lines
5 poke55,.:poke56,56:clr
10 dv=peek(186):ifdv<8thendv=8
15 poke53280,0:poke53281,0:print"[147]"
16 poke53371,0
25 poke53272,31
30 ad=49152
35 sysad:sysad+12
100 rem rigorous precession
105 rem (j2000.0)
110 p1=(NULL):r1=p1/180
112 bx$="[158]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_[157]_"
113 sb$="[153]_[157]_[157]_[157]_[157]_[157]_[157]_"
114 tq=ti
115 e0=2000:a$="a"
120 print"[147]":sysad+9,6:poke53272,31
121 print"[153] ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
122 printtab(3)sb$;""tab(33)sb$
123 print"[153] ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
124 poke214,2:print:printtab(9)"[129][211][212][193][210] [204][207][195][193][212][207][210] [208][210][207][199][210][193][205]!"
125 poke214,4:print:printtab(5)"[158][201]nitial [197]poch [155](yr) ";:l9%=4:gosub670:ie=q9
130 print:printtab(5)"[158][198]inal [197]poch [155](yr):";:l9%=7:gosub670:fe=q9
135 ifabs(ie-fe)<=500 then150
140 poke214,10:print:printtab(9)"[212]ime span too long!"
141 tq=ti
142 ifti<tq+200then142
145 goto120
150 print"[147]":sysad+9,17:gosub415
151 print"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
152 printbx$;""tab(38)bx$
153 print"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^"
155 ifa$="n"then185
160 print"":printtab(9)"[153] [201]nitial [197]poch [146][158]";ie
165 printtab(3)"[154][210].[193]. ";
170 gosub350:a=v
175 printtab(3)"[150][196]ec ";
180 gosub350:d=v
185 ma=0:md=0:gosub570
190 a2=a1:d2=d1
195 ifa$="n"then220
200 printtab(9)"[153] [208]roper [205]otion in "
205 printtab(3)"[159][210].[193]. [155](sec/yr) ";:l9%=7:gosub670:m1=q9
210 printtab(3)"[159][196]ec. [155](arcsec/yr) ";:l9%=7:gosub670:m2=q9
215 printtab(1)"[158]^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^":sysad+9,8
220 ma=m1:md=m2:gosub570
225 v=a1:gosub380
230 printtab(9)"[153] [198]inal [197]poch [146][158] ";fe
235 printtab(3)"[159][210].[193]. [155](h,m,s):[158] ";
240 printv1;v2;v3
245 v=d1:gosub380
250 printtab(3)"[159][196]ec. [155](d,m,s):[158] ";s$;
255 printv1;v2;v3
260 a3=(a1-a2)*3600/ny
265 d3=(d1-d2)*3600/ny
270 printtab(9)"[153] [208]roper [205]otion in "
275 printtab(3)"[159][210].[193]. [155](sec/yr):[158] ";
280 printa3
285 printtab(3)"[159][196]ec. [155](arcsec/yr):[158] ";
290 printd3
295 print
300 printtab(9)"[159] [211]elect [207]ne [208]lease:[146] "
305 printtab(3)"[150]([193])[153]nother [211]tar [150]([206])[153]ew [198]inal [197]poch"
315 printtab(11)"[150]([213])[153]niverse [205]enu
320 [161]a$:[139]a$[179][177]"a"[175]a$[179][177]"n"[175]a$[179][177]"u"[167]320
325 [139]a$[178]"a"[167]150
330 [139]a$[178]"n"[167]110
335 [139]a$[178]"u"[167]3050
350 [143] input sexagesimal
355 s[178]1:[153]"conthour/defdegree (def+/-):";:l9%[178]3:[141]670:q9[178][197](q9$):v$[178]q9$
356 [153][163]9)"contminutes:";:l9%[178]4:[141]670:v2[178]q9
357 [153][163]9)"contseconds:";:l9%[178]3:[141]670:v3[178]q9
360 [139][200](v$,1)[178]"-"[167]s[178][171]1
365 v1[178][182]([197](v$))
370 v[178]s[172](v1[170]v2[173]60[170]v3[173]3600)
375 [142]
380 [143] output sexagesimal
385 s$[178]"+":[139]v[179]0[167]s$[178]"-"
390 v[178][182](v):v1[178][181](v)
395 vm[178]60[172](v[171]v1):v2[178][181](vm)
400 v3[178]60[172](vm[171]v2)
405 [142]
410 :
415 [143] pression parameters
420 ny[178]fe[171]ie
425 t0[178](ie[171]e0)[173]100:t1[178]ny[173]100
430 t2[178]t1[172]t1:t3[178]t1[172]t1[172]t1
435 h1[178]2306.2181:h2[178]1.39656
440 h3[178].000139:h4[178].30188
445 h5[178].000345:h6[178].017998
450 k1[178]1.09468:k2[178].000066
455 k3[178].018203
460 l1[178]2004.3109:l2[178][171].8533
465 l3[178][171].000217:l4[178][171].42665
470 l5[178][171].000217:l6[178][171].041833
475 w[178](h1[170]h2[172]t0[170]h3[172]t0[172]t0)[172]t1
480 zt[178]w[170](h4[170]h5[172]t0)[172]t2[170]h6[172]t3
485 zd[178]w[170](k1[170]k2[172]t0)[172]t2[170]k3[172]t3
490 th[178](l1[170]l2[172]t0[170]l3[172]t0[172]t0)[172]t1
495 th[178]th[170](l4[170]l5[172]t0)[172]t2[170]l6[172]t3
500 zt[178]zt[172]r1[173]3600:zd[178]zd[172]r1[173]3600
505 th[178]th[172]r1[173]3600
510 [143] zt,zd,th = euler amgles
515 :
520 [143] rotation matrix
525 s1[178][191](zt):c1[178][190](zt)
530 s2[178][191](zd):c2[178][190](zd)
535 s3[178][191](th):c3[178][190](th)
540 xx[178]c1[172]c3[172]c2[171]s1[172]s2
545 yx[178][171]s1[172]c3[172]c2[171]c1[172]s2:zx[178][171]s3[172]c2
550 xy[178]c1[172]c3[172]s2[170]s1[172]c2
555 yy[178][171]s1[172]c3[172]s2[170]c1[172]c2:zy[178][171]s3[172]s2
560 xz[178]c1[172]s3:yz[178]s1[172]s3:zz[178]c3
565 [142]
570 [143] proper motion correction
575 a0[178](a[170]ma[172]ny[173]3600)[172]15[172]r1
580 d0[178](d[170]md[172]ny[173]3600)[172]r1
585 :
590 [143] sherical rectangular
595 sa[178][191](a0):ca[178][190](a0)
600 sd[178][191](d0):cd[178][190](d0)
605 x0[178]ca[172]cd:y0[178]sa[172]cd:z0[178]sd
610 [143] 3-d transformation
615 x1[178]x0[172]xx[170]y0[172]yx[170]z0[172]zx
620 y1[178]x0[172]xy[170]y0[172]yy[170]z0[172]zy
625 z1[178]x0[172]xz[170]y0[172]yz[170]z0[172]zz
630 [143] rectangular sperical
635 a1[178][193](y1[173]x1)
640 [139]x1[179]0[167]a1[178]a1[170]p1
645 [139]a1[179]0[167]a1[178]a1[170]2[172]p1
650 a1[178]a1[173](r1[172]15):[143] final r.a.
655 d1[178][193](z1[173][186](x1[172]x1[170]y1[172]y1))
660 d1[178]d1[173]r1:[143] final dec.
665 [142]
670 q9$[178]"":[151]198,.
675 [161]a$
680 [151]646,[187](1)[172]15[170]1:[153]"*cmd";:[139]a$[178]""[167]675
685 [139]a$[178][199](13)[167][153]" ":q9[178][197](q9$):[142]
690 [139](a$[178][199](20)[175][195](q9$))[167]q9$[178][200](q9$,[195](q9$)[171]1):[137]720
695 [139][195](q9$)[177][178]l9%[167][158]ad[170]9,1:[137]675
700 [139](a$[177][178]"0"[175]a$[179][178]"9")[176]a$[178]"."[176]a$[178]"-"[176]a$[178]"+"[167]710
705 [137]675
710 q9$[178]q9$[170]a$
715 [153]"sys"a$;:[158]ad[170]9,7:[137]675
720 [153]" cmdcmd cmd";:[137]675
3050 [158]ad[170]15
3060 [153]"loadstopload"[199](34)"b.universe iii"[199](34)","dv
3070 [153]"run28"
3080 [151]631,13:[151]632,13:[151]198,2:[128]
10000 d[178][194](186):n$[178]"0:star locator":[159]15,d,15,"s"[170]n$:[160]15:[148]n$,d:[128]